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ORIGINAL ARTICLE 

A Semi-Mechanistic IVIetabolism IVIodel of CYP3A 
Substrates in Pregnancy: Predicting Changes in 
IVIidazolam and Nifedipine Pharmacokinetics 

SK Quinney' ^ AN Mohamed^ MF Hebert", DWl Haas^ ^ S Clark«, JG Umans'^ SN Caritis' and L Li^°; for the Obstetric-fetal 
Pharmacology Research Unit Network 

Physiological changes in pregnancy, including changes in body composition and metabolic enzyme activity, can alter drug 
pharmacokinetics. A semi-mechanistic metabolism model was developed to describe the pharmacokinetics of two cytochrome 
P450 3A (CYP3A) substrates, midazolam and nifedipine, in obstetrics patients. The model parameters were optimized to fit the 
data of oral midazolam pharmacokinetics in pregnant women, by increasing CYP3A-induced hepatic metabolism 1 .6-fold in the 
model with no change in gut wall metabolism. Fetal metabolism had a negligible effect on maternal plasma drug concentrations. 
Validation of the model was performed by applying changes in volume of distribution and metabolism, consistent with those 
observed for midazolam, to the pharmacokinetics parameters of immediate-release nifedipine in healthy volunteers. The predicted 
steady-state areas under the concentration-time curve (AUCs) for nifedipine were within 1 5% of the data observed in pregnant 
women undergoing treatment for preterm labor. This model predicts the pharmacokinetics of two CYP3A substrates in pregnancy, 
and may be applicable to other CYP3A substrates as well. 
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The results of database surveys in developed countries show 
that 44-93% of pregnant women take at least one prescription 
medication other than vitamins and iron J Pregnancy leads to 
a number of physiologic changes that alter drug metabolism 
and disposition. Despite the impact of pregnancy on pharma- 
cokinetics, this population is underrepresented in and often 
excluded from clinical trials. Therefore, for many drugs, the 
pharmacokinetics and effects in pregnant women are largely 
unknown. One means of overcoming the ethical and practical 
constraints associated with studying the effects of individual 
drugs in pregnant women is to use pharmacometric models 
to predict drug disposition at various stages of pregnancy. 

Traditionally physiology-based pharmacokinetic (PBPK) mod- 
els incorporate physiologic parameters (such as blood flow and 
tissue volumes) with drug-specific parameters (such as physi- 
cochemical properties and in vitro metabolism data) to predict 
the concentration of drugs in various tissues. Such PBPK mod- 
els of pregnancy have historically been used by toxicologists to 
predict the concentration of toxicants, including drugs, in various 
tissues."'^ More recently pharmacokinetic modeling has come 
to be accepted by regulatory agencies and pharmaceutical 
scientists as an important component of clinical pharmacology 
studies that address drug efficacy and drug-drug interactions." 



We have proposed a semi-mechanistic pharmacokinetic 
model of drug metabolism, based on standard one- or two- 
compartment models and incorporating both hepatic and 
intestinal metabolism.'"^ This approach uses compartmental 
pharmacokinetic parameters available from the literature or 
derived from clinical study data. Metabolic enzyme activity 
can be estimated from in vitro rate constants or from clinical 
data. This model has been validated for the cytochrome P450 
3A (CYP3A) probe substrate, midazolam, using data from 
drug-drug interaction studies with ketoconazole, clarithromy- 
cin, diltiazem, and erythromycin in healthy volunteers.'"^ 

Here, we adapt this semi-mechanistic model to describe 
drug disposition in pregnancy, incorporating placental and 
fetal compartments and physiologic changes associated with 
pregnancy (Figure 1 and Table 1). The parameters of the 
model are optimized using midazolam clinical pharmacoki- 
netic data in pregnancy;^" sensitivity analyses are performed 
to illustrate the impact of fetal metabolism, changes in hepatic 
blood flow, plasma protein binding, and CYP3A activity on 
plasma midazolam concentrations. Based on the adjusted 
parameter estimates determined for midazolam, the model 
was used for predicting the concentration levels of another 
CYP3A substrate, nifedipine, in pregnant women. 
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Figure 1 Obstetric semi-mechanistic metabolism model. Central and 
peripheral volumes of distribution are derived from empirical two- 
compartment pharmacokinetics analyses. Clearances in gut wall, liver, 
and fetal liver (CL^j^,, CL^, CL,^,^,) are estimated by the well-stirred model. 
CLp, renal clearance; k^, absorption rate constant; Op,,, portal vein blood 
flow; 0„, hepatic blood flow; 0^^, hepatic artery blood flow; 0^,^^ placental 
blood flow; Q,^,, fetal (umbilical vein) blood flow. 



Table 1 Model parameter values" 





Midazolam, 
healthy" 


Midazolam, 
pregnant' 


Nifedipine, 
healthy" 


Nifedipine, 
pregnant" 


Vc(l) 


43 (28%) 


45.9 (28%) 


20(28%) 


23(28%) 


Vp (1) 


88.4 (20%) 


93.9 (20%) 


30 (20%) 


35.5 (20%) 




1.17(40%) 


2.5 (40%) 


3.36(40%) 


4.48(40%) 


f 


0.04 


0.048 


0.05 


0.06 




2,755 (30%) 


4,592(30%) 


1 ,245 (30%) 


1 ,992 (30%) 




2,549 (30%) 


2,549 (30%) 


1 ,245 (30%) 


1 ,245 (30%) 


(\imo\n) 


5.8 


5.8 


5.0 


5.0 






0.65 




0.65 


Ke,u.(l) 




3.4 




3.4 






88 




167 


(^mol/l) 




114 




34 






10.8 




10.8 






15.9 




15.9 



All other parameters as in Quinney etal." and provided in supplementary 

information online. 

■Data presented as mean (%CV). 

RESULTS 

Postpartum midazolam model 

The healthy volunteer model for midazolam^"^ was used to 
simulate plasma midazolam concentrations after a 2-mg 
oral dose of the drug. Predicted plasma midazolam con- 
centrations were similar to those observed in women 6-10 
weeks postpartum (Figure 2a). ^° Observed and predicted 
areas under the concentration-time curve (AUCs) from time 
0 to oo(AUC^) were 17.4 and 16.9 ng-h/ml, respectively 
(residual sum of squares, 20.8; percentage error, 3%). 

IVIidazolam obstetric model 

The initial semi-mechanistic obstetrics model for midazolam 
metabolism was based on the model derived from data in 



healthy human volunteers, with the addition of placental and 
fetal compartments (Figure 1), a 20% increase in fraction 
unbound (f ),^° a twofold increase in intestinal and hepatic 
CYP3A clearance ( V^gJ," and an increased volume of distri- 
bution {V^Y^ based on reports in the literature. In comparison 
with observed midazolam data at 28-32 weeks' gestation,^" 
the model underestimated (AUCq) (5.7 vs. 9.5 ng-h/ml) and 
maximal plasma concentration (C^^^) (2.1 vs. 6.4ng/ml) 
(Figure 2b and Table 2). 

Sensitivity analyses 

In order to identify the impact of individual parameters on 
the midazolam concentration-time curve, select parameter 
values were varied 0.5- and 2-fold relative to the initial model 
(Figure 3 and Table 2). Changes in hepatic blood flow (O^) 
and the absorption rate constant (/cj had little effect on AUG, 
but a large effect on C^^^. A twofold change in had large 
effects on the plasma concentration-time curve. Similarly, 
reducing V or V ^ increased both the AUG and the 

^ max.GW max.H 

C of midazolam. 

max 

The contribution of the fetus to maternal drug pharmacoki- 
netics was also examined using sensitivity analysis. Fetal 
hepatic V^^^ was estimated to be 88mg/h, based on in vitro 
GYP3A7 metabolism of midazolam and estimated expres- 
sion of GYP3A7 in the fetal liver.^^" Gompletely excluding 
fetal clearance (V = 0) did not alter the maternal con- 

\ max,fet ' 

centration-time course of midazolam. Even increasing the 
Knaxfet '^V much as 20-fold reduced midazolam AUG in 
maternal plasma only minimally (Table 1). Although fetal 
metabolism has little effect on maternal plasma drug con- 
centration, it may impact fetal exposure. However, given the 
low intrinsic clearance of midazolam by the fetal liver, even 
a complete exclusion of fetal metabolism is predicted to 
increase fetal (AUG^) of midazolam by only 2%. Addition- 
ally, changes in fetal blood flow (Q^^ and 0^,^^ did not sig- 
nificantly affect plasma concentrations (see Supplementary 
Figure S1 online). It therefore appears likely that the effect 
of fetal metabolism on maternal plasma midazolam concen- 
tration is minimal. 

The volume of distribution of midazolam is expected to be 
higher in pregnancy because of increased plasma volume, 
increased extracellular volume, and decreased binding. We 
assumed that in pregnancy, the increase in the central volume 
of distribution is proportional to the change in plasma volume, 
and that peripheral volume of distribution accounts for addi- 
tional increases in the total volume of distribution. Sensitivity 
analysis indicates that the maternal plasma drug concentra- 
tion-time profile was not sensitive to changes in the volume 
of the peripheral compartment (see Supplementary Figure 
S2 online). 

The original model, which anticipated that the effects of 
pregnancy on GYP3A V^^^ would be similar in the gut wall 
and the liver, led to underestimation of C^^^. Based on results 
of the sensitivity analysis and iterative model fitting, the initial 
model of midazolam disposition in pregnancy was revised 
such that /c was increased to 2.5 h-\ V ^ was increased 

a ' max.H 

1 .6-fold, and V^^^^^ was left unchanged from the healthy vol- 
unteer model. This modified pregnancy model was able to 
accurately predict midazolam AUG within 1%, but continued 
to underestimate C^^ by 23% (Figure 2c). 
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Figure 2 Predicted and observed midazolam concentrations, (a) Predictions of postpartum concentrations of midazolam by the hiealtfiy volunteer 
semi-mechianistic metabolism model, (b) Initial predictions by the obstetric model based on empiric changes assumed from the literature, (c) Visual 
predictive check of final obstetric midazolam model. Observed data are from women at 26-30 weeks of gestation. '"The participants received a single 
2 mg dose of midazolam orally. Median and 5th and 95th prediction intervals are indicated by black and gray lines, respectively. 



Table 2 Effect of varying model parameters on AUG, C^^^^, and model fit (residual sum of squares) 

Residual sum of 





Fold-change" 


AUC (ng h/ml) 


% Error for AUC 


'-max (na/ml) 


% Error for C^^^ 


squares 


Initial model" 




5.7 


-40 


2.1 


-67 


31.2 


Q„(l/h) 


0.5x 


5.7 


-40 


1.6 


-75 


44.1 




2x 


5.8 


-39 


2.5 


-60 


22.6 




0.5x 


2.8 


-70 


1.2 


-81 


17.1 




2x 


11.5 


20 


3.4 


-46 


52.4 




0.5x 


5.7 


-40 


1.4 


-78 


51.7 




2x 


5.7 


-40 


3.0 


-53 


15.8 


K„.x,ow(mg/h) 


0.5x 


8.4 


-12 


3.1 


-51 


15.4 




2x 


3.5 


-63 


1.3 


-79 


51.2 




0.5x 


11.4 


20 


3.4 


^6 


17.1 




2x 


2.8 


-70 


1.2 


-81 


52.4 


K„.x,«(mg/h) 


0 


5.7 


^0 


2.1 


-67 


31.2 




20x 


5.5 


-42 


2.1 


-67 


31.2 


Final model" 




9.5 


-0.9 


4.8 


-23 


4.3 


Observed" 




9.5 




6.4 






AUG, area under the concentration-time curve; 
"Fold-change in parameter from initial model. "Q^ = 
f^ = 0.048, k^ = 2.5 h-i ; V^^^ =2,549 mg/h, V^^„ = 


^, maximal plasma concentration. 

°92.3 l/h, f =0.048, ;c^ = 1.17h"'; l/,„a„=5,09 

4,592 mg/h, l/,^^ = 88mg/h.''Hebert ef a/." 


8 mg/h, V„^,„ = 


5,510mg/h, l/^^^,^, = 88mg/h.=0„ 


= 92.3 l/h. 



Application to nifedipine 

In order to validate the increased liepatic CYP3A activity 
associated with pregnancy, we applied our model to another 
CYP3A substrate, nifedipine. Using data from the literature, 
we constructed a model of nifedipine pharmacokinetics 
in healthy volunteers (Table l).'^''*' The two-compartment 
parameters and in vitro metabolism kinetics adequately 
predicted the pharmacokinetics of oral nifedipine (10mg) in 
healthy volunteers (Figure 4a). The mean predicted AUC(,_gOf 
nifedipine, 156ng-h/ml (95% confidence interval of 149-163), 
was similar to the mean observed AUC(,_g, 180ng-h/ml (95% 
confidence interval of 1 50-21 0). Parameter adjustments were 
made as in the midazolam model (Table 1), and the steady- 
state concentration of nifedipine after oral administration of 
a 20 mg immediate-release capsule every 6h was predicted 
in pregnant subjects (Figure 4b). As compared to observed 
data in subjects taking nifedipine for preterm labor,'^ the 
median AUG was underpredicted by 11% (210 (121-299) 



ng-h/ml vs. 237 (224-253) ng-h/ml), and C^^^ was underpre- 
dicted by 3% (178 (166-188) ng/ml vs. 184790-308) ng/ml). 

DISCUSSION 

Pharmacokinetic studies of the CYP3A substrates, midazo- 
lam and nifedipine, have recently been conducted in pregnant 
women (25-34 weeks of gestation).'"^' Both drugs are predom- 
inately metabolized by CYP3A, the intrinsic clearance (CLJ 
showing CYP3A4 > CYP3A5 > CYP3A7.i^ Unlike many other 
CYP3A substrates, neither drug is significantly transported 
by p-glycoprotein.^**''^ Clinical studies^"" indicate that CYP3A 
activity increases in pregnancy. By incorporating various physi- 
ologic changes of pregnancy into a pharmacokinetic model, we 
attempted to isolate the effect of pregnancy on CYP3A activity. 

PBPKs models have long been used to estimate fetal toxici- 
ty'' However, relatively few models have been developed to esti- 
mate the pharmacokinetics of therapeutic drugs in pregnancy. 
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Figure 3 Sensitivity analysis of obstetric semi-mectianistic metabolism model. The solid line indicates predicted maternal plasma drug concentration 
based on the initial obstetric semi-mechanistic metabolism model. Broken and dotted lines represent twofold decreases and increases, respectively, 
in (a) hepatic blood flow, Q^; (b) absorption rate constant, k^; (c) fraction unbound, f ; (d) intestinal rate of metabolism, V^^^^; and (e) hepatic rate 
of metabolism, V^^^^^. (f) Effect of fetal metabolism on maternal plasma concentrations as predicted by the initial model (solid line), assuming no 
contribution from fetal metabolism (dotted line), and with a 20-fold increase in cytochrome P450 3A7 (CYP3A7) V^^ (dashed line). 



The primary purpose of this model is to estimate changes in 
drug clearance and maternal pharmacokinetics of drugs. We 
have simplified the model to focus on important organs of drug 
metabolism, namely, Intestine and liver. Although not impact- 
ing the pharmacokinetics of these CYP3A substrates, placen- 
tal aromatase and (perhaps) fetal hepatic enzymes may play 
a role In the metabolism of other drugs (e.g., methadone as a 
substrate of CYP19). Therefore, we believe that separation of 
fetal and placental compartments may be relevant In the case 
of some drugs. This Is similar to models proposed by Gabriels- 
son and Paalzow.^ We have additionally collapsed the PBPK 
model Into two compartments: peripheral distribution and cen- 
tral distribution. This is In contrast to PBPK models proposed 
by Gaohua etal. and Andrews etal., which propose 13 and 20 
maternal compartments, respectlvely.^"'^^ While our model may 
not fully describe pregnancy-related changes In drug distribu- 
tion, this simplification allows us to focus on the roles of the 
relevant organs during this Initial stage of model development. 

Model development 

The original semi-mechanistic drug metabolism model was 
developed based on midazolam pharmacokinetics In healthy 
volunteers taking part in drug-drug interaction studies.'"" The 
healthy volunteer model Incorporated the key organs involved 
In the absorption, metabolism, and distribution of oral drugs: 
Intestinal lumen, Intestinal wall, liver, portal circulation, and 
systemic circulation. Other physiologic compartments, which 
are not key determinants of drug clearance, were collapsed 
Into the central and peripheral compartments. Compart- 
ments associated with drug metabolism are assumed to be 
part of the fast-equlllbrating compartment, and their volumes 
are therefore deducted from that of the central compartment. 



Drug metabolism is described seml-mechanlstlcally, based 
on the well-stirred model, enabling the interrogation of fac- 
tors that may affect drug pharmacokinetics (e.g., blood flow, 
drug binding, and enzymatic activity). The original model was 
validated by its adequate prediction of plasma concentration- 
time curves for midazolam in postpartum settings (Figure 2a). 
It also indicates that, by 6-10 weeks postpartum, the factors 
controlling midazolam pharmacokinetics have largely returned 
to those expected In the absence of pregnancy. 

In order to describe the pharmacokinetics of midazolam 
In pregnancy, the original model was altered to Incorporate 
physiologic changes associated with pregnancy: additional 
placental and fetal compartments with metabolic potential, 
Increased volume of distribution, decreased plasma protein 
binding, and Increased CYP3A clearance. The Initial model 
assumed that both Intestinal and hepatic GYP3A activity 
Increased twofold." However, It underestimated observed 
^max (Figure 2b). Several factors can affect C^^^, Including 
bioavailability, absorption rate, volume of distribution, and 
drug clearance. Changes In systemic clearance may result 
from changes in hepatic blood flow, protein binding, or CYP3A 
Intrinsic clearance. Using sensitivity analyses and knowledge 
of physiologic changes In pregnancy, we explored the Impact 
of each of these factors on midazolam disposition. 

Sensitivity analysis 

Intestinal absorption of drugs may be affected in pregnancy 
by alterations In gastric motility and pH.^^ Being a hlgh-per- 
meabllity drug, midazolam is nearly completely absorbed 
from the gastrointestinal tract. However, changes In motility 
can impact the rate of intestinal absorption (I.e., kj, altering 
both C and time to reach C (t ). Sensitivity analyses 

max max ^ max' ^ ^ 
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Figure 4 Model validation with nifedipine, (a) Plasma concentration-time 
curve for nifedipine in liealthy volunteers after a single lOmg oral dose, 
(b) Steady-state plasma concentration-time profile of nifedipine 10mg 
given orally every 6 h to women at 24-34 weeks of gestation. Solid lines 
represent the 5th and 95th percentiles of the simulations (n - 1,000). 
The broken line represents the median of the simulations. Open and 
closed circles represent the observed nifedipine plasma concentrations 
of (a) nine studies in healthy volunteers, and (b) individual observations 
in pregnant women, respectively. 

indicated that increasing tine of midazolam improved pre- 
diction by the obstetric model (Figure 3b). 

Tliere are some data, tliough limited, to suggest that preg- 
nancy lias little effect on hepatic blood flow.^^ In addition, the 
sensitivity analysis results suggested that changes in hepatic 
blood flow would have little effect on the AUG of oral midazo- 
lam (Figure 3a). Therefore, it does not appear that changes 
in hepatic blood flow are responsible for the increased clear- 
ance of midazolam during pregnancy. 

Midazolam is highly bound to albumin, and therefore even 
small changes in protein binding are capable of producing 
marked changes in the drug's systemic clearance. This is 
shown in the sensitivity analysis, where a decrease in f from 
0.048 to 0.024 has an effect on midazolam pharmacol<inetics 
similar to that resulting from a twofold increase in hepatic V^^^ 
(Figure 3 and Table 2). In pregnancy, albumin concentra- 
tions decrease by -20%, with average albumin concentra- 
tions in the third trimester being 38.8 ± 2.9 g/l as compared 
to 47.7 ± 2.5 g/l in nonpregnant women. Consequently, the 
protein binding of drugs that bind primarily to albumin, such 
as midazolam and nifedipine, is expected to decrease, lead- 
ing to increased clearance. However, the effect of protein 



binding cannot fully account for the twofold increase in clear- 
ance observed in pregnant women. 

Given that changes in liver blood flow and protein binding 
do not completely account for increased midazolam clear- 
ance, it is possible that altered metabolism, either in the 
liver or fetal/placental unit, may play a role. The presence 
of CYP3A activity in the placenta is still a matter of debate. 
CYP3A4, CYP3A5, and CYP3A7 mRNA have been identified 
in preterm placenta, but the amounts are much lower than 
those found in the liver.^""^' However, only one study^' has 
observed CYP3A activity in placenta. It is unknown whether 
other placental enzymes (e.g., CYP19) are capable of mida- 
zolam metabolism. Because of the incompleteness of current 
knowledge, we did not consider placenta as a potential driver 
of midazolam clearance in our model. Fetal liver expresses 
CYP3A7 to a greater extent than CYP3A4 or CYPSAS.^" 
In vitro studies of midazolam with CYP3A isoforms,^^ and 
knowledge of CYP3A7 expression in the fetal liver, provided 
the data for estimating the contribution of the fetal liver to 
midazolam disposition. As shown in Figure 3f, metabolism of 
midazolam in the fetal liver had a negligible effect on mater- 
nal plasma midazolam concentrations. 

Consequently, it appears that the increased midazo- 
lam clearance in pregnancy is largely caused by increased 
hepatic CYP3A activity. Interestingly, even a twofold increase 
in CYP3A activity in the model, both in the gut wall and in 
the liver, led to an underestimation of C and AUC. Assum- 

' max 

ing a 20% increase in f^, hepatic CYP3A activity appears to 
increase 1.6-fold in pregnancy. Independent regulation of 
CYP3A activity in the liver and intestines is consistent with 
results of studies in mice, which showed that hepatic CYP3a 
activity increased 2.2- and 2.5-fold on gestational days 15 
and 19, respectively, but that no change was observed in 
intestinal CYP3a activity or expression.^" However, studies in 
the nonhuman primate Macaca nemestrina™ have not found 
increases in hepatic or intestinal CYP3A expression or mida- 
zolam clearance; the mechanism underlying the increase in 
hepatic CYP3A activity in humans is unknown. 

Iterative analyses indicated that maintaining V^^^^ at non- 
pregnancy values while increasing V^^^^, f , and from the 
initial healthy volunteer parameter estimates resulted in the 
best fit of the observed pregnancy data. In order to achieve a 
better fit of the absorption profile, k^ was estimated to increase 
to 2.5h-\ As this increase in absorption rate was unexpected, 
we further interrogated the postpartum data. Changing k^ in the 
healthy volunteer model from 1.17 to 2.5 h"^ appeared equally 
capable of predicting the postpartum data (see Supplementary 
Figure S3 online). Therefore the increase in k^ appears to be 
attributable to the sample population, and is not an effect of 
pregnancy. 

Validation with nifedipine 

The pregnancy-induced changes proposed by the obstetric 
semi-mechanistic metabolism model were validated using 
nifedipine. Two-compartment pharmacokinetic parameters for 
nifedipine administered to healthy volunteers were obtained 
from the literature^^'^** and adapted to fit the semi-mechanistic 
metabolism model by incorporation of in vitro and V^^ (see 
Supplementary Data online). Predictive checl<s of the model 
against published concentration-time profiles (Figure 4) 
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indicated tliat the liealthy volunteer model described the data 
adequately. Changes in the volume of distribution, protein 
binding and hepatic CYP3A activity equivalent to those in the 
obstetrical midazolam model were applied. The model predic- 
tions of the mean nifedipine concentrations in subjects taking 
20 mg nifedipine for preterm labor were within 1 5% of observed 
data. Several concentration values near the early points were 
overpredicted by the model, possibly because of the high inter- 
individual variation observed in absorption profiles.^^ 

Obstetric PBPK models of pregnancy have typically 
focused on toxicokinetics and risk assessment.'' However, 
physiologic and drug-specific parameters can be incorpo- 
rated into pharmacokinetics models to enhance understand- 
ing of drug metabolism in pregnancy. The eventual goal of 
such work is to develop a "virtual pregnant woman" model 
to predict the pharmacokinetics and effects of drugs in preg- 
nancy. However, such models require the incorporation of 
gestational-age-dependent physiologic changes. While some 
of these physiologic changes can be estimated on the basis 
of gestational age or weight-based algorithms, quantita- 
tive knowledge of other changes is limited and more diffi- 
cult to acquire. Further in vitro and in vivo studies of drug 
metabolism and transport during pregnancy are required to 
fully populate these virtual models. Future development of 
this model will incorporate gestational-age-related growth 
changes. Additionally, the effects of pregnancy on drug distri- 
bution and absorption, as well as other routes of elimination, 
will be incorporated into future models. 

A limitation of this and other pharmacokinetic models is the 
difficulty in model validation. In practice, drug concentrations 
can be determined only in the central compartment (plasma). 
Pharmacokinetic models enable prediction of drug concen- 
trations in other compartments, e.g., at the site of drug action 
or in the fetus. However, validation of these concentrations 
is unlikely to be carried out in humans. Fetal concentration 
values can be obtained noninvasively at birth from umbili- 
cal cord blood. However, these values are limited to a single 
time point after maternal dosing, and may not reflect drug 
exposures at earlier ages of gestation. To complicate matters, 
both mother and fetus are exposed to drugs taken by the 
mother, requiring consideration of both therapeutic effective- 
ness and the potential for toxicity. 

In conclusion, we have shown that our midazolam semi- 
mechanistic metabolism model predicts midazolam pharma- 
cokinetics in women in a postpartum setting. By incorporating 
physiologic changes associated with pregnancy, we show 
that increasing midazolam hepatic clearance 1 .6-fold in the 
model leads to accurate estimates of the clearance of oral 
midazolam in pregnant women. This was confirmed using 
data from a second CYP3A substrate, nifedipine. However, 
further study would be required to fully validate the lack of 
change in intestinal clearance despite increased hepatic 
clearance. For example, a pharmacokinetics study involving 
semi-simultaneous intravenous and oral administration of 
a probe drug would be informative. The application of this 
model to substrates metabolized by other enzymes will allow 
us to better understand the changes in drug pharmacokinet- 
ics that are associated with pregnancy. This model is a first 
step in the development of predictive models for individual- 
ized dosing in pregnancy 



METHODS 

Semi-meclianistic metabolism midazolam pliarmacoidnetic 
model. The semi-mechanistic metabolism pharmacokinetic 
model for midazolam in nonpregnant healthy subjects has 
been reported previously^-'' For the obstetric model, two 
additional compartments were added to represent the pla- 
centa and fetus (Figure 1). 

Elimination of midazolam occurs in the maternal gastro- 
intestinal tract, liver, and kidney. In addition, the placenta 
and fetal liver have drug metabolizing capabilities. However, 
placental activity of CYP3A isoforms is minimal, ^''-'^^ and no 
evidence exists to indicate that midazolam is metabolized by 
the placenta. Therefore the placenta was not considered an 
organ of elimination in our model. 

Maternal hepatic metabolism occurs through CYP3A-medi- 
ated and other enzymatic clearance pathways, e.g., glucuroni- 
dation. The CYP3A- and non-CYP3A-mediated intrinsic 
clearances of midazolam in the liver are estimated as: 



ci.„, 



Vriax,3A 



Vriax, nonSA 



(1) 



(2) 



where V ,„ and V 

max,3A max,non3A 



are the maximum velocity of 
metabolism by CYP3A and non-CYP3A mechanisms, is 
the IVIichaelis-IVIenten constant, f is the unbound fraction of 
midazolam, and C„ is the hepatic drug concentration. 

Intrinsic clearance of midazolam within the gut wall (CL^j^,,^) 
is assumed to be entirely by CYP3A; therefore. 



ci.„, 



VTiax,3A 



(3) 



where V is the maximum velocity of metabolism 
by CYP3A in the gut wall. The value for the gut wall is 
assumed to be equivalent to hepatic K^. 

The fetal liver primarily expresses CYP3A7. The metabo- 
lism of midazolam in the fetal liver, as in the adult liver, can 
be estimated by scaling up in vitro V^^^ and values to esti- 
mate V , using Equation (4):'''' 



K„ax,,e, = K„ax,3A7 xlSEFxCYP3A7 abundanco 



xMPPGLx Liver weight 



(4) 



where V is the maximal reaction velocity in the fetal liver. 



V 



is the maximal velocity determined in rCYP3A7 in 



vitro, ISEF is an intersystem extrapolation factor,'"' CYP3A7 
abundance is the amount of CYP3A7 enzyme per mg 
microsomal protein, and IVIPPGL is the amount of microsomal 
protein in mg per gram of liver. Williams et al. reported the 
in vitro V^^ and values of midazolam hydroxylation by 
recombinant CYP3A7 as 4 nmol/min/nmol and 114 jjmol/l, 
respectively.^" It is estimated that fetal liver microsomes con- 
tain 0.3 nmol CYP450 protein per mg.^" Based on ICRP data, 
the fetal liver weighs 130g at 38 weeks gestation. MPPGL is 
estimated to be 26mg/g.'"' In the absence of data, ISEF was 
assumed to be 1 . Therefore, 1/^^^^^, is estimated as 88mg/h. 
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Volumes of distribution established from compartmental 
pliarmacol<inetic analyses are used to estimate the volumes 
of the central and peripheral compartments. Liver, portal vein, 
gut, placenta, and fetal volumes were based on reported 
physiologic volume estimates, and these volumes are 
subtracted from the central distribution compartment. During 
pregnancy, women gain 12.5l<g of body weight on average. 
In the semi-mechanistic metabolism midazolam model, the 
portion of weight gain not accounted for by the fetus and pla- 
centa is divided between the central and peripheral compart- 
ments. The assumption is that extracellular fluid and blood 
volume reside in the central compartment, and that changes 
in other tissues and fat stores are a component of the periph- 
eral compartment. 

While total cardiac output increases in pregnancy, the 
value of is assumed to be similar to that in the nonpreg- 
nant state, and is calculated as an allometric expression of 
total body weight:^^ 

Qh = 3.75 X (body weight)"-"' (5) 

Hepatic artery flow (O^J and portal venous blood flow 
(Opy) represent 25 and 75% of O^, respectively Placental 
(*^piac) ^'"^ umbilical venous (QJ blood flows were estimated 
at 30 weeks of gestational age.^** 

Calibration data. Hebert et ai. determined maternal plasma 
midazolam concentrations after a single 2mg oral dose 
in women at 28-32 weeks of gestation.^" The study was 
repeated in the women at least 6 weeks postpartum. Indi- 
vidual plasma concentration-time profiles of these subjects 
were compared with the predicted concentration-time pro- 
files for 500 individuals. 

Sensitivity analysis. Sensitivity analysis was conducted to 
examine the effects of V^^^,, V^^^^, Q,, V^, and fetal 
metabolism {V^^^^J on the model. Starting from the ini- 
tial model, each parameter was individually increased and 
decreased twofold, except for V , , which was reduced to 0 

' ~ max,fei 

and increased 20-fold (Table 2). 

Parameter optimization. Based on results of the sensitiv- 
ity analysis, individual parameters or subsets of them were 
manually adjusted without taking interindividual variability into 
account, and holding all other parameters constant. Predicted 
concentrations were compared with the mean concentration- 
time curve data reported by Hebert ef a/.^°The optimal model 
was chosen based on visual inspection of the concentration 
vs. time plot and minimization of the residual sum of squares. 

Nifedipine semi-mechanistic metabolism model. A healthy 
volunteer model and an obstetric semi-mechanistic metabo- 
lism pharmacokinetic model were constructed for another 
CYP3A substrate, nifedipine, using the same model struc- 
ture as for midazolam. Central and peripheral volumes of 
distribution and were determined from published two-com- 
partment analyses of nifedipine pharmacokinetics in healthy 
volunteers (Table 1).^'-^^ Nifedipine plasma concentrations in 
healthy volunteers after oral administration of 10mg imme- 
diate-release capsules were extracted from mean concen- 
tration values reported in eight different studies involving 



healthy volunteers, using Grab-It Graph Digitizer (Datatrend 
Software, Raleigh, NG).^'^ In vitro V^^^ and values of nife- 
dipine in human liver microsomes were obtained from pub- 
lished studies (see Supplementary Data online). All other 
physiologic parameters, i.e., blood flows and physiologic 
organ volumes, were identical to the midazolam healthy vol- 
unteer semi-mechanistic metabolism model. 

Nifedipine obstetric semi-mechanistic metabolism model. The 
pregnancy-associated changes in CYP3A activity, volume, 
and plasma protein binding identified in the midazolam semi- 
mechanistic metabolism obstetric pharmacokinetic model 
were applied to the corresponding model for nifedipine. Given 
that protein binding of midazolam and nifedipine are similar, 
changes in f and central and peripheral volumes of distribu- 
tion were applied as per the final midazolam model. Hepatic 
V was increased 1 .6-fold to account for increased CYP3A 

max 

activity in pregnancy. 

Nifedipine validation data. Predicted concentration-time 
curves of nifedipine after 20 mg oral dosing to pregnant women 
were compared to data obtained as part of the Under-studied 
Drugs in Pregnancy protocol conducted by the Obstetric- 
Fetal Pharmacology Research Unit network.^" The study 
examined nifedipine pharmacokinetics in 20 women (24-36 
weeks gestation) who were receiving 10-20mg immediate- 
release nifedipine every 4-8 h for preterm labor, as clinically 
indicated. Given that the majority of women were taking 20 mg 
every 6 h, the 1 0 mg doses were normalized to 20 mg for visu- 
alization and comparison. The study was approved by the 
institutional review board at each study site (Georgetown Uni- 
versity, University of Pittsburgh, University of Texas Medical 
Branch-Galveston, and University of Washington), and written 
informed consent was obtained from all the participants. 

Software and model analysis. All simulations were performed 
in R V 2.13.1,"'' using the odesolve package. Noncompart- 
mental analysis was performed using the PK package"" 
to obtain AUC(,_^ (midazolam) or AUGo_j, (nifedipine). The 
example code is provided as online supplementary informa- 
tion (see Supplementary Data online). Predictions were 
compared with observed data using visual predictive checks 
and calculation of percentage error ([predicted - observed/ 
observed] and residual sum of squares. 
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Study Highlights 

WHAT ISTHE CURRENT KNOWLEDGE ON THE 
TOPIC? 

Physiologic changes in pregnancy, including changes 
in body composition and metabolic enzyme activity, can 
greatly affect drug pharmacokinetics, and thereby effec- 
tiveness and safety of drugs. Physiology-based pharma- 
cokinetic models have traditionally been applied in fetal 
toxicology studies. However, such models are also capa- 
ble of predicting therapeutic activity. 

WHAT QUESTION DIDTHIS STUDY ADDRESS? 

The aim of this study was to develop a semi-mechanistic 
metabolism model for CYP3A substrates in pregnancy. 
This model can help identify mechanistic changes 
responsible for the observed increased clearance of 
CYP3A substrates in pregnancy. 

WHATTHIS STUDY ADDSTO OUR KNOWLEDGE 

Using our model for CYP3A substrates in pregnancy, 
we are able to identify potential mechanisms for the 
pregnancy-associated changes in midazolam and nife- 
dipine pharmacokinetics. Assuming a 20% increase 
in f^, the model indicates that hepatic, but not intes- 
tinal, metabolism by CYP3A increases 1.6-fold during 
pregnancy. 

HOWTHIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY ANDTHERAPEUTICS 

This study supports the use of pharmacometric modeling 
to identify mechanisms of pharmacokinetic changes in 
pregnancy and provides a framework for future obstetric 
pharmacokinetic models. 
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